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Abstract 

This talk is a review of our studies of instantons and their properties as 
seen in our lattice simulations of SU (2) gauge theory. We have measured 
the topological susceptibility and the size distribution of instantons in the 
QCD vacuum. We have also investigated the properties of quarks moving 
in instanton background field configurations, where the sizes and locations 
of the instantons are taken from simulations of the full gauge theory. By 
themselves, these multi-instanton configurations do not confine quarks, 
but they induce chiral symmetry breaking. 
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"Nonperturbative QCD-Structure of the QCD Vacuum" 



1 Introduction 



What features of the QCD vacuum are responsible for confinement or for the 
generation of the observed structure of hadron spectroscopy? This question 
might, in principle, be answered by lattice simulations of non-Abelian gauge 
theories. We have been studying the properties of instantons from SU (2) lattice 
simulations. This talk is a survey of our work as presented in Refs. jl], g| ||, |J, || . 

The study of the QCD vacuum using lattice Monte Carlo is complicated by 
two problems. The first one is that the dominant features of the QCD vacuum 
as seen in lattice simulations are short distance fluctuations (as they would be 
for any quantum field theory). They are basically uninteresting noise. The so- 
lution to this problem is to invent operators which are insensitive to the short 
distance behavior of the field variables. This brings the second problem: The 
separation of vacuum structure into short distance and long distance parts is 
ambiguous, and what one sees can depend strongly on the operators one uses. 
All direct smoothing transformations^) distort the original lattice configuration. 
This makes the extraction of (continuum) short to medium distance physics, like 
observations of topological objects, very delicate. If any space-time symmetric 
smoothing transformation is repeated enough times, all the vacuum structure 
in any finite volume, including the simulation volume, will be washed away. 
Thus, it does not make sense to extrapolate one's results to the limit of a very 
large number of smoothing steps. The only measurements which are physically 
meaningful are measurements which are extrapolated back to the original lat- 
tice, that is, back to zero smoothing steps. This requires careful monitoring of 
observables over the whole history of smoothing transformations. 

In QCD instantons may be responsible for breaking axial symmetry and 
resolving the U{1) problem. The relevant observable is the topological sus- 
ceptibility \t, defined as the infinite volume limit of 



where Q is the topological charge and V the space time volume. In QCD Xt 
is a dimension-4 object with no weak coupling expansion, and a calculation of 
Xt in physical units in the continuum requires nonperturbative techniques. In 
the large- N c limit the mass of the rj' is (probably!) related to the topological 



The left hand side of this equation is equal to (180 MeV) 4 in the real world. 
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Based on phenomenological models, it has been argued that instantons are 
largely responsible for chiral symmetry breaking and the low energy hadron and 
glueball spectrum. || [T^] Instanton liquid models attempt to reproduce the 
topological content of the QCD vacuum and conclude that hadronic correlators 
in the instanton liquid show all the important properties of the corresponding 
full QCD correlators. These models appear to capture the essence of the QCD 
vacuum, but their derivations involve a number of uncontrolled approximations 
and phenomenological parameters. 

Instanton physics on the lattice is as full of controversy as continuum instan- 
ton physics. There are presently three different ways of measuring a topological 
charge. The "geometric" definition |Q, [T^| reconstructs a fiber bundle from the 
lattice gauge field and identifies the second Chern number of this bundle with 
the topological charge. It will always give an integer, but if the configuration 
is sufficiently rough, it can fail catastrophically. "Algebraic" definitions [^3) in- 
troduce some lattice discretization of Q, as a sum of closed paths of loops. The 
worst aspect of the algebraic definitions is that the topological charge can mix 
with quantum fluctuations. Finally, one can define Q through fermionic opera- 
tors. (For a recent example, see Ref. All these definitions have a cutoff: 
they cannot see instanton (like) configurations when the instanton radius is too 
small (typically p/a ~ 1 — 2). It is usually not possible to specify this cutoff 
precisely, and it can contaminate lattice measurements. 



2 Finding Instantons 

We have explored three different methods for extracting information about 
topology from lattice simulations. Each has its own strengths and weaknesses. 



2.1 Inverse Blocking 

"Inverse blocking" does not distort the original configuration. This technique 
has been introduced by Hasenfratz|l5| in his lectures. Imagine beginning with 
a set of lattice variables {V} on a lattice whose spacing is a and lattice size is 
L. The lattice action is a fixed-point (FP) action |H| S FP (U). The inverse 
blocking transformation constructs a set of fine lattice variables {U} occupying 
a lattice of lattice spacing a/2 and lattice size 2L, by solving the steepest-descent 
equation 

S FP (V) = min (S FP (U) + kT(U, V)) , (3) 
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where kT(U, V) is the blocking kernel. Inverse blocking identifies the smoothest 
among the configurations that block back to the original configuration. Since 
for fixed point actions, topology is unchanged by inverse blocking, the strategy 
is to take a rough configuration (generated by Monte Carlo using a FP action), 
inverse block it, and measure the charge on the fine lattice, where the instantons 
are twice as big and charge measurement is more reliable. We did thisQ and 
found Xt — (235(10) MeV) 4 for SU(2). This created a low-level controversy, 
since our number was larger than other measurements. [l8], |l9f 

Inverse blocking is reliable, but it is very expensive. In four dimensions, 
one can only perform it once (to go from lattice spacing a to a/2). The fine 
configurations arc still too rough to identify individual instantons. 



2.2 Cycling 

Our smoothing mechanism for seeing instantons is called "cycling." [|| One 
first performs an inverse blocking from a coarse lattice to a set of fine lattice 
variables by solving Eqn. |[ Now the original lattice occupies one of the 16 
sublattices of the fine lattice. Next, we perform a blocking transformation to a 
set of coarse variables {W} based on one of the other sublattices. The delicate 
coherence among the fine variables is broken and the new coarse variables are 
strongly ordered on the shortest distance scale (as measured, for example, by 
the expectation value of the plaquette) while retaining all long distance physics 
(because they are generated by a RG blocking transformation). This is the 
second part of the cycling transformation V^(x) — » Uu{x) — * W fl (x). Cycling 
steps can be iterated, and a few cycling steps can reduce the plaquette to within 
0.001 of its free-field value. 

Individual instantons can be seen after a few cycling steps. Their sizes drift 
with smoothing (see the next section for pictures) but the drift is small enough 
that we can extrapolate their properties back to zero cycling steps. We found 
Xt^ = 230(10) MeV for SU(2), an instanton density of about two per fm 4 , and 
a mean instanton radius of about 0.2 fm. Again, the susceptibility is higher than 
others', and now the mean instanton size is smaller than other measurements, 
although it agrees with the predictions of the interacting instanton liquid model. 




Cycling is still quite expensive because of the inverse blocking step. This 
makes it hard to push to small lattice spacing and to test scaling. 
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2.3 RG Mapping 



Cycling produces a sequence of coarse ({V}) and fine ({£/}) lattices {V} 
{Ui} — > {Vi} — > {C/2} — * {V2} — > RG mapping is a technique for eliminating 
the inverse blocking step and generating a sequence of coarse lattices {V} — ► 
{Wi} — > {VF 2 } — > . . . where {W n } is an approximation to {V n }. The idea is 
that, while formally the inverse blocking is non-local, for local FP actions the 
dependence of the fine links on the original coarse links dies away exponentially 
with their separation Q and the mapping {V} — > {Ux(V)} — > {Vi(V)} can be 
considered local. The {W„} lattice can be constructed from the original coarse 
{V} lattice as a sum of loops of {V}'s which are designed to reproduce a cycled 
{V n }. We discovered that this could be done by APE-smearing: Q from a set 
{V} construct a new set of links {X} by 

X„(x) = (1 - c)V„{x) + c/Q^{V v {x)V^x + v)V^{x + *>) f 

+ V 1/ (x-P)%(x-i>)V li (x-P + (i)), (4) 

with X^x) projected back onto SU(2) to generate W^(x). We found that a 
series of steps with c = 0.45 was a good choice to mock up cycling. 

RG-mapping is very cheap, but the price is that everything which is mea- 
sured must be monitored carefully, and extrapolated (if necessary) back to zero 
mapping steps. 

The problem is that any approximation to the fixed point charge is expected 
to distort the charge density profile. We have to correct this distortion. This 
can be done by monitoring how the charge density-measured by the FP charge- 
changes in the course of smoothing and extrapolating this back to zero smooth- 
ing steps. Of course we cannot directly measure the FP charge since this would 
involve several inverse blocking steps. However after a few smoothing steps the 
configurations become smooth enough so that our improved charge operator is 
very close to the exact FP charge and data taken after further smoothing can 
be used for extrapolating back. 

As an example, we consider a 16 4 configuration generated with Wilson action 
j3 = 2.5. Fig. [l] shows the size of the 8 "stable" objects as a function of the 
APE-smearing steps. From the 4 instantons (diamonds) 3 increase in size while 
one decreases, but all vary linearly with smearing steps. The slope of the linear 
change for all of them is small. Three of the anti-instantons (bursts) behave 
similarly, though one has a slightly larger slope. The fourth anti-instanton 
(crosses) starts to grow rapidly after 18 smearing steps and will disappear after 
a few more steps. This object is likely to be a vacuum fluctuation, not an 
instanton. 
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Figure 1: Radius versus APE-smearing steps of instantons (diamonds) and 
anti-instantons (bursts and crosses) on a 16 4 (3 — 2.5 configuration. 



tn the early stages of cycling (1-2 steps) we see many lumps of charge. 
Most of them quickly grow or shrink away. The locations of what we call true 
topological objects are stable over many smearing steps and their size changes 
slowly. To identify them on the lattice one has to track them over several 
smearing steps and monitor their behavior. 

Instantons present in QCD simulations differ from hand crafted instantons 
in trivial background configurations in that the former usually grow while the 
latter objects always shrink under APE smearing. 

Neither cycling nor RG-mapping affect the string tension, but the short 
distance part of the potential is distorted. 

One of our goals was to compare results obtained with the cycling/RG map- 
ping method with results published using other algorithms (l7|, [l^, pj| , so we 
used the Wilson action in conjunction with RG-mapping. 

We observe, as expected, a small systematic decrease in the susceptibility 
as we increase the number of smearing steps. At large [3 the change is small 
and statistically insignificant. Only at j3 — 2.4, where the configurations are the 
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roughest, do the 12 and 24 smearing steps results differ by about a standard 
deviation. 

The susceptibility increases by 10% from (3 = 2.4 to (3 = 2.5 but stabilizes 
after that at the value \t ^ = 220(6) MeV. We interpret the change between the 
extrapolated (3 = 2.4 result and the larger (3 results as due to the absence of small 
instantons at (3 = 2.4 because of the larger lattice spacing. This interpretation 
will be supported by the instanton size distribution result discussed below. 

We have identified individual instantons after every 2 APE-smearing steps 
between 12 and 24 steps. Since these configurations are still rough, many of the 
objects identified as instantons are in fact vacuum fluctuations and disappear 
after more smoothing steps. Figure |^ shows the observed instanton size distri- 
bution on the 12 and 24 times smeared lattices at f3 = 2.5. For comparison 
we also plot the result of Ref. jlTj which corresponds, in our normalization, to 
about 100 times smeared lattices. It is obvious from the figure that the total 
number of identified objects decreases as we increase the smoothing. The den- 
sity of identified objects is 4.6 per fm 4 after 12 smearing steps, 3.0 per fm 4 after 
24 smearing steps, and about 2 per fm 4 in Ref. fl7| . These density values are 
considerably larger than the expected value of about 1 per fm 4 . The maximum 
of all 3 distributions is around p « 0.3 fm, but that does not mean that on 
the original lattice p ~ 0.3 fm since instantons usually grow under smearing. 
This growth can also be observed from the increasing tail of the distribution, 
especially after the many blocking steps of Ref. [jljj . 

Our final result for the instanton size distribution (after extrapolation back 
to zero smoothing steps) is shown in Figure [| where we overlay the data ob- 
tained at [3 = 2.4 (octagons on 12 4 , squares on 20 4 lattices), (3 = 2.5 (diamonds), 
and (3 = 2.6 (bursts). Since the smoothing method cannot identify instantons 
with p < 1.5a, we chose the bins such that the second bin for each distribution 
starts at p = 1.6a. That means that we expect the second bin of each distribu- 
tion to be universal. The first bins on the other hand contain only some of the 
small instantons and their value is not expected to be universal. 

The four distributions form a universal curve indicating scaling. The (3 = 2.4 
curves cover only the p > 0.2 fm region, and small instantons are obviously miss- 
ing. The agreement between the 12 4 and 20 4 configurations at (3 — 2.4 indicate 
that a linear size of about 1.4 fm is sufficient to observe all the topological 
objects. The (3 = 2.5 and 2.6 distributions have most of the physically rele- 
vant instantons, supporting the scaling behavior observed for the topological 
susceptibility 

The instanton liquid model predicts a very similar picture to ours. In 
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Figure 2: The size distribution on 12 (bursts) and 24 (crosses) times APE 
smeared lattices at /3 = 2.5. The square symbols are the results of Rcf. 11 
rescaled appropriately. 
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Figure 3: The size distribution of instantons. Octagons correspond to (3 = 2.4 
12 4 , squares to j3 = 2.4 20 4 , diamonds to (3 = 2.5 and bursts to (3 = 2.6. The 
first bin of each distribution is contaminated by the cut-off. The solid curve is 
a two parameter fit to the data points according to the formula in Ref. 23. The 
dashed curve is a similar fit from Ref. 23 which describes the instanton liquid 
model quite closely. 
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Ref. [M Shuryak predicted an instanton size distribution that peaked around 
p = 0.2 fm. The curve in Fig. || is a fit using the two loop perturbative instanton 
distribution formula with a "regularized" log 

8ir 2 

S I = ——=b a L + b l \ogL (5) 
9 2 {P) 

L = -log[(pA inst )- p + C p ] (6) 

where bo and b\ are related to the first two coefficients of the perturbative (3 
function and p and C are arbitrary parameters. The solid curve is a fit to 
our data while the dashed curve is the fit given in Ref. |H| which describes 
the instanton liquid model quite closely. The difference between the two fits is 
significant at large p values only. We do not know if changing the parameters 
of the interacting instanton liquid model slightly would change the predictions 
of that model improving the agreement with the Monte Carlo data. 

To make a long story short, we believe that there are two reasons why our 
instanton numbers are different from those of Refs. jl7], |lj, First, we are 
sensitive to instantons of smaller radius. As we cut larger and larger instantons 
from our sample, we see Xt f&h- Second (and most important) we extrapolate 
our measurements back to zero smoothing steps. Instantons generally grow 
under smoothing, so our extrapolated sizes are smaller. 



3 What Do Instantons Do? 



To test what role instantons play in the QCD vacuum, we took a set of SU(2) 
configurations with lattice spacing a ~ 0.14 fm and smoothed them enough 
to identify the instantons. These smoothed configurations have essentially the 
same string tension as the original configurations, even though about seventy 
per cent of their vacuum action is carried by the instantons. We then identified 
the sizes and locations of the instantons in the configurations and built multi- 
instanton configurations from them. We built both "parallel" and randomly 
oriented instanton configurations. 

Notice that the instanton locations are not random, and the sizes are not 
taken from a model distribution-they come from the (lattice) QCD vacuum. 

In Fig. [| the heavy-quark potentials obtained from the three ensembles are 
compared. We can conclude that neither the parallel nor the randomly oriented 
instantons confine. It seems that instantons by themselves are not responsible 
for confinement. 
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Figure 4: The heavy-quark potential measured on the 9 times cycled real config- 
urations (octogons), the randomly (crosses) and the parallel (squares) oriented 
instantons. 
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Figure 5: (am^) 2 vs. mo the bare quark mass for the staggered action, with 
raw gauge links (diamonds) and 9-cycled links (squares). The lightest mass in 
the pseudoscalar channel with instanton background configurations is shown by 
octagons. 



Next, we computed hadron spectroscopy in these background configurations. 
One could think of spectroscopy on the smoothed configurations as spectroscopy 
computed with a new kind of quark action, one which is insensitive to the short 
distance behavior of the gauge field. This is similar in spirit to the use of "fat 
links" by the MILC collaboration |22| and by J.-F. Lagae and D. K. Sinclair 
[ p3| , or of the approximate FP actions developed by one of us.p4[ 

Staggered fermions provide the most interesting results. It appears that 
cycling does not affect the p mass vs the Goldstone pion. However, on the 
smoothed lattices the non-Goldstone partner of the pion is degenerate with the 
pion within observational uncertainty. On the original lattices the two mesons 
have a mass ratio of about 1.4. 

The purpose of this exercise is not to test spectroscopy in detail. The im- 
portant point is that the action of the many-times-cycled configurations is dom- 
inated by instantons, and yet the asymptotic heavy quark potential and hadron 
spectroscopy are basically unchanged from what we saw on untouched configu- 
rations. 
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Figure 6: The pseudoscalar propagator from 9-cycled configurations, with stag- 
gered fermions of bare mass am® = 0.05. 

But these results do not answer the question: Are instantons responsible 
for long distance physics, or are the structures responsible for long distance 
physics some other objects carrying low action, which have been preserved by 
the smoothing transformation? 

We computed spectroscopy on the two types of multi-instanton ensembles 
using both staggered and Wilson fermions. The dominant feature of both spec- 
troscopy calculation is that the quarks are deconfined. This is seen most easily 
in the staggered fermion spectroscopy. 

Consider the pseudoscalar propagator of the 9-cycled configurations, shown 
for one quark mass in Fig. ^|. It looks like any generic lattice pseudoscalar, a 
more-or-less pure hyperbolic cosine with no oscillations. The staggered fermion 
pseudoscalar propagators on instanton background fields are quite different: 
they show the characteristic sawtooth pattern of free antipcriodic fermions. 
(See Fig. 0). 

Fig. H shows the lattice for staggered fermions on the 9-cycled and in 

instanton background configurations, (ipip) m the instanton background tracks 
the value of {?pip} measured on the 9-cycled configurations quite closely, down to 



13 




2 Q I I I I I I I I I I I I I I I I 

5 10 15 

t 

Figure 7: The pseudoscalar propagator in (randomly rotated) instanton back- 
ground configurations, with staggered fermions of bare mass amo — 0.05. The 
curve is a fit to a single propagating particle plus the qq branch cut. 
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Figure 8: (tpt^) from raw configurations (diamonds), 9-cycled configurations 
(squares) , and instanton-background configurations which are parallel (crosses) 
and randomly oriented (octagons), vs. bare quark mass amy. The line shows 
the free-field value, 2mo- 

small quark mass. It appears that the instantons, present in equilibrium gauge 
field configurations of the QCD vacuum generated by Monte Carlo, are breaking 
chiral symmetry by themselves. This effect is a cornerstone of instanton-liquid 
models of hadron structure. 

If (tpip) is nonzero, one expects the spectrum contains a would-be Goldstone 
boson (the pion) in addition to massive quarks and (possibly) other resonances. 
To test this hypothesis, we fit the pseudoscalar correlator to two terms: a pure 
hyperbolic cosine (a pole in the frequency plane), plus a qq branch cut, with 
the quark mass as a parameter. We used the analytic expression for the branch 
cut (expressed as a momentum and frequency mode sum) , with the appropriate 
boundary conditions and source used in the simulations. In the random instan- 
ton background, we clearly see a light mass which decreases towards zero as the 
quark mass vanishes. (See Fig. ||.) 

The quark mass in the randomly oriented instanton background is also de- 
termined by the fit. It is close to the bare mass and is less than half the rho 
mass-and less than half the pion mass- (as measured on the smoothed lattices). 
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Are there other bound states? Only for the randomly rotated instanton 
configurations could we see a convincing pseudoscalar, and so we restricted our 
analysis to that data set. We again attempted to fit to a single resonance plus 
a qq continuum. Of course, this assumption is questionable. The fits are not 
of high quality and our results should not be taken too seriously: In the "SC" 
channel (saturated by the 7r 2 and scalar mesons in the confined phase) we saw 
a light bound state whose mass roughly tracked the mass of the (presumed) 
pion resonance in the pseudoscalar channel. (That is, there is a multiplct of 
Goldstonc bosons.) In the vector channels, a state with a mass 750-850 MeV 
appears, in addition to the free qq continuum. This state has about the same 
mass as the vector meson in the confined system. However, the dominant feature 
of all these channels is still the free qq continuum, with fitted quark masses of 
100 MeV or lower whose energy is always lower than the mass of any (presumed) 
resonance. 

The reader should be aware that the lattice spacing is large enough that it 
could contaminate the coupling of quarks to instantons. 



4 Conclusions 

We have measured the instanton content of the SU(2) lattice vacuum. In- 
stantons do not confine, but they seem to be connected to chiral symmetry 
breaking. We are presently using RG-mapping to study the properties of in- 
stantons in SU(3), both in the pure gauge theory and for full QCD. This spring 
we expect to begin computing hadron spectroscopy in smoothed configurations 
and in multi-instanton background configurations. How similar will they be? 
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